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Abstract 

In the context of the Loop Quantum Cosmology we have analysed the holonomy correction to 
the classical evolution of the simplified Bianchi I model in the presence of vector fields. For the 
qq Universe dominated by a massive vector field or by a combination of a scalar field and a vector 

{S| field a smooth transition between Kasner-like and Kasner-unlike solutions for a Bianchi I model 

has been demonstrated. In this case a lack of initial curvature singularity and a finite maximal 
energy density appear already at the level of General Relativity, which simulates a classical Big 
13^ Bounce. 

Introduction 



Vector fields play a prominent role in modern particle physics. We detect vector fields in laboratory 
and we expect them to participate actively in the physics of the early Universe. As shown in 
literature, vector fields can contribute to generation of the observed large scale structure [U El El HI [5] , 
CSj background anisotropics [BJ, inflation [7j and gravitational waves [SJ, so that, in principle, one could 
^vq reconstruct standard FRW cosmology without scalar fields. In order to obtain slow-roll evolution and 
flat power spectra of initial inhomogeneities, vector fields need to have a time dependent mass term. 
It may originate from the Higgs mechanism, non-canonical kinetic terms or from a non-minimal 
• • coupling to gravity. 

On the other hand, a study of the influence of Loop Quantum Cosmology correction on evolution 
X of FRW universes Efl QH Q21 EE] and Bianchi I models O HQ QH HSJ HH] has been performed. In 
both cases the existence of the quantum of length changes significantly the evolution of scale factors 
near the Planck scale, resolving the problem of the initial singularity. 

The Big Bounce model usually assumes, that the Universe is filled with a homogeneous scalar 
field, so there is no physical argument for breaking the 0(3) symmetry of the metric tensor. In this 
paper we shall analyse the LQC in the Universe filled with a vector field. We start with a simple 
example, where the Universe is dominated by a massless or a massive vector field. Next, to study 
the influence of the vector field in more realistic set-ups we consider a case, where a vector and a 
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scalar fields are present. They can be uncoupled, or strongly coupled. The last case corresponds in 
this paper to the time dependent mass term for the vector field. An interesting idea is to introduce 
vectors with a non-minimal coupling to gravity, however it is not known how to take into account 
the LQC corrections in this case. 

In all of those scenarios the vector field produces anisotropic pressure, which forms a physical 
source for the anisotropic expansion in the Bianchi I model. The main purpose of this paper is 
to find out, whether vectors could change significantly the behaviour of scale factors close to the 
moment of the Big Bounce. 

In this paper we investigate two generalisations of the isotropic LQC model. Two schemes of 
constructing the pi functions are presented and compared. It is worth to note, that both of them 
give similar results. 

This paper is organized as follows: In section [T] we introduce the Hamiltonian formalism for 
the metric tensor with two scale factors and we calculate the Einstein's equation together with 
equations of motion for matter fields. In section [2] we analyse the influence of the loop correction on 
the classical evolution of scale factors and matter fields for the case of the vector field domination. 
In section [3] we study the Universe dominated by both, vector and scalar fields in the context of 
LQC. Conclusions are presented in section [4] 



1 The classical theory 

1.1 Ashtekar formalism and equations of motion 

In this section we introduce Ashtekar variables for the anisotropic Universe. The full discussion 
of this issue can be found in |14| 120] . Let us consider a homogeneous, massive vector field pQ 
= A(t)S 3 ^, which contributes to the Universe's energy density. Since the vector field breaks the 
0(3) symmetry of the Universe, one needs to introduce two scale factors - b(t) for the direction 
pointed by the vector and a(t) for transverse directions. This means, that the model is a simplified 
version of diagonal Bianchi I. Let us set the lapse function to be N(t) = 1, which means, that the 
metric tensor is given by 

ds 2 = dt 2 - a 2 {t){dx 2 + dy 2 ) - b 2 {t)dz 2 = Diag(l, -a 2 , -a 2 , -b 2 ) . (1.1) 



Using the eq. (1.1) one can introduce the canonically conjugated Ashtekar variables [H 

pi = L\L 2 ab , p 2 = L\a 2 , c\ = Lijd , c 2 = L 2 r yb , (1.2) 

where 7 is the Barbero-Immirzi parameter and {c\,pi} = 87r 2 7 , {c 2 ,p 2 } = 8irG"f. L\,L 2 are 
comoving lengths introduced to make Ashtekar variables independent of the rescaling of scale factors. 
One constructs the fiducial cell V Q = L\L 2 , which is used to obtain finite scalar constraint. Ashtekar 
variables pi have the following interpretation: the area perpendicular to the z axis is proportional to 
a 2 and areas perpendicular to x an y axes are proportional to ab. On the other hand variables c\ and 
c 2 have the interpretation of external curvature along x, y and z axes respectively. Together with 
Ashtekar variables we shall also define the total physical volume of the fiducial cell V = {p\p2) 1 ^ 2 ■ 



For the metric tensor described by (1.1) one obtains V oc a b, so V is proportional to the volume of 
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the Universe. Let us assume, that the Universe is filled with a vector field and a scalar fielcQ With 
the convention 8irG = 1 the scalar constraint for the minimally coupled general relativity looks as 
follows 

1 3//2 1 

C = — 2 i={c\p\ + 2cic 2 ;pip 2 ) + ^vri + ^-m 2 A 2 + — tt^ + pxjpiV {<!>) = , (1.3) 

TPly/pl 2p 2 ' 2 Pi 2piy^i v 



where 7ta = P^J 2 A/px = o?A/b and ir^ = p\^/p~2 are canonical momenta of vector and scalar fields. 
The scalar constraint decomposes into C = C grav +C m at, where C ma t = Px^fp~2P- Equations of motion 
for the Ashtekar variables are given by Hamilton equations n 

■ / rx 7 dC ■ s rx dC n A^ 

PX = {PX,C} = --— , P2 = {P2,C} = -J— , (1.4) 
2 OCi OC2 

ci = {ci,C} = ~~ , c 2 = {c 2 ,C} = 7^- , (1.5) 

2 Opi Op2 



which together with the eq. (1.3) give Einstein's equations and the equations of motion for the 
matter fieldfE 

H (H + 2%) = p = P(f> + PA , (1.6) 

- + ^ + HU = -p ± = - P4> - p A , (1.7) 

2^ + H 2 = -p {{ =- P(j> +pA , (1.8) 
P(j> + {2H + H){ P ^+p 4> ) = , Pa + 2H(pa+pa)+'H(pa-Pa) = , (1.9) 

where H = - ,H = |. These equations, together with equations of motion for matter fields, give 
the complete evolution of the classical space-time. In realistic cosmological models one expects 
the vector field to produce small background anisotropies of the Universe, so one could expect 
A = H — % <C H + %. Nonetheless, as we shall see later the vector field domination leads to 
observable anisotropies. 

Let us focus for a while on the background evolution of the vector field. One can show, that 

PA = ^2 + ™ 2A % PA = 7^(A 2 - m 2 A 2 ). (1.10) 

The pressure and the energy density are physical quantities, so they are invariant under the re- 
scaling of any scale factor. This means, that A transforms like A — > l z A under the rescaling 
b — > l z b. Then A(t) shall be replaced by the physical field U(t) = A/b, which gives ^A 2 = U 2 and 
^A 2 = U 2 + 2UUU + U 2 U 2 . The eq. Esty looks now as follows 

U + (2H + n)U + (2-HH + n + m 2 )U = 0. (1.11) 



In stead of a scalar field one can consider a perfect fluid. 
2 Coefficients 7 and 7/2 come from the symmetry reduction of diagonal Bianchi I model. They appear because 
{ci,pi} = 7/2 and {c 2 ,p 2 } = 7. 

3 We refer to scalar and vector fields collectively as matter fields. 
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One can see, that U has time-dependent effective mass term m 2 = 2HH + 1-1 + m 2 . For a = b one 
obtains in 2 = 2H 2 + H + m 2 = R/6 + rn 2 , so U shall be very heavy in the high energy regime, in 
which H 2 m 2 . During inflation the subdominant U is much heavier than the inflaton, so after 
a few e-folds it reaches its minimum and starts oscillating. The energy density of a vector field 
behaves matter-like, pa oc a" 3 oc e~ 3Ht , and vectors are expected to be washed out during inflation. 
On the other hand the power spectrum of the vector field perturbation would not be flat enough to 
fit the CMB data. For the time of the cosmic inflation we need thus an extra tachyonic mass term 
m 2 ~ —2H 2 which would allow the slow roll evolution. After the reheating the vector field could 
play a role of the curvaton and dominate the Universe via its oscillations. One can learn more about 
the issue of the vector field effective mass from e.g. [H El [7]. 

One shall remember, that the Ashtekar formalism gives us the same results as the Lagrange 
formalism with the action 



S 



d X\ 



-9 



l R - l F uv F^ + l m 2 A^A, 
2 4 M 2 ' 



where R is the Ricci scalar of the anisotropic Universe 

R 



n d b, 

2 {H 2 + 2HU + 2- + -) 
a b 



:i.i2) 



Cl-13) 



1.2 Massless vector field domination 

The most popular case in the LQC is a scalar field without any potential term. This produces simple 
solution and allows a matter field to play a role of a time variable. Since it is the most fundamental 
calculation in LQC we will first of all consider this case for the vector field in the GR frame. As 
we will see later on this model will bring us rather non physical solutions of Einstein equations. 
On the other hand those solutions shall describe the low energy limit for the massless vector field 
in the LQC frame. When the energy-stress tensor consist only of the massless vector field, then 
PA = P = P- From eq. ( 1.6|1.8 ) one obtains 



H 1 + 2HU 



2- + fr 
a 



Ed 



J1.14) 



a b 
a b 

where E = const > has a dimension of time. This shows, that the Universe shrinks along the z 



axe if only a < 0. On the other hand, the eq. (1.9) gives 

p + 4Hp = => p = pa 



a 



J1.15) 



where p a = p(t Q ) and t is any fixed moment in the Universes history. This equation can be also 
obtained from the equation of motion of the vector field. For the massless vector field one gets 



A + (2H - H)A = => A = exp 



(H - 2H)dt 



Then, from the eq. (1.14) one obtains 



A 



a 



D 2 al R 2 
2a 2 b 2 



2^2 



D z a 
2^ 



:i.i6) 



11.17) 
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combine eq. ( 1.8|1.15 ) which together give 
This equation can be simplified to 



where D = const > 0. From the eq. (1.15) one obtains D = ^/2p a E. To calculate a(t) let us 

(1.18) 



2da + a 2 = p at a 2 



a 2 

Po 2 



:i.i9) 



where C\ = const has a dimension of energy. One can see, that this solution faces a constraint. Let 
us define the initial value of the scale factor a by aj = a(tj) = a C\j p Q . Then a(ai) = b(ai) = 0, 
which means, that to obtain real values of a, b one needs a > aj. Thus the energy density cannot be 
greater than pj = p (a / 'a/) 4 . It is important to note, that the initial (maximal) energy density is 
simply a constant of integration. So far let us make no assumptions about pj, but later on we shall 
investigate a massless vector in the LQC frame, for which the natural scale of the maximal energy 
density is the Planck scale. Thus, pi will set by hand to that scale. Formulae presented so far take 
the simplest form for t = tj, which is the choice we shall adopt. 

So far we have calculated H,%,b } p,A as functions of the scale factor a. Thus, to find time 



dependence of those quantities one needs to find a solution of the eq. (1.18), which is 

3 



a 



+ 3 



a 



aj 
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w 



pi(t-t T y 



[1.20) 



By shifting the time coordinate t —> t + ti one obtains w 



pit. The exact solution of (1.20) is 



a(t) = ai 



2 i/3 



(2 + w 2 + V^w 2 + w 4 )V3 



+ 



(2 + w 2 + V4w 2 + w 4 ) 1/3 
21/3 



:i.2i) 



From eq. ( 1.13|l.l4|l.21 ) one obtains R = for any t. For t — > one also obtains R^ U R^ U — > Apj 



and R^ VOL ^ R^vafi —> 20p|. Thus the massless vector field domination gives no initial curvature 
singularity of the Univers e j 4 ^} which together with aj ^ 0, H(tj) = are features of the Big-Bounce. 
On the other hand from ( 1 . 14|1 .20 ) one obtains % — > oo, b — > for t — > 0, which looks alike to the 
Big Bang scenario. When we are very close to the initial value of the energy density one obtains 
w<l, which gives 



ai 



1 + 



ai 



i + - Pl r 



b oc t 



V oc t 



i + - Pl r 



n.22) 



The Universe expands along the z direction for w E (0,4). Let us note, that a{t) for w<Cl evolves 
like a(t) around the Big Bounce in the dust dominated isotropic Universe in LQC frame For 
p <^L pi one obtains w S> 1. Then 



a 



2/3 



ai 



1/3 



t 2/3 



b oc t 



-1/3 



oc a 



-1/2 



V oc t 



[1.23) 



4 The energy conditions for the vector field domination require p > 0, p > and p > \p\. The massless vector field 
satisfies them, since p = p — A 2 /2b 2 > 0. 

5 Let us note, that the Big Bounce in the GR frame may be also obtained in i.e. cycling Universe [2T] and f{R) 
gravity [22]. 
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The low energy limit bring us to the Kasner-like solution of Bianchi I model, with one direction 
shrinking and two expanding. The model with such predictions is obviously excluded by the as- 
tronomical observations, but as we have mentioned it is just a GR limit of the Universe with the 
LQC correction, which shall be discussed later. The evolution of a(w) and b(w) is shown in the fig. w 



1.3 Massive vector field 

In this subsection we shall investigate the case of a massive vector field domination with m = const. 
In our calculations in GR frame we wanted to construct the low-energy limit for the LQC solutions. 
Since then we are interested in the massive vector field with purely kinetic initial conditions similar 



to the massless field scenario 1.2 If initially Aj = then we can describe the beginning of a Universe 



by the massless vector field domination. Thus, from the eq. (1.17) one obtains 



a i\ 1 /?. — it. 2 tt /?, — a i f-, a i 



A = JlpiEai 1 - — ~ -JlpiEaipiF, V = J2p 1 —[\- — \ (1.24) 
V a / 4 a V a J 

After a short period of this analytically described evolution we need to take into account correc- 
tions form the mass term. Then, for the numerical simulation we shall assume a finite %(i*)and 
a non-zero H(t+), where * denotes the boundary condition for analytical and initial condition for 
numerical calculation. Regardless of any initial conditions one obtains % > for e.g. t > O(10 3 i p j) 
and m = 10 -4 . When the energy density density decreases, and m > max{H,7i}, the vector field 
starts to oscillate and one obtains the matter-like dominated Universe. The becomes pressure- 
less, so there are no anisotropies in the energy-stress tensor. Then, as it is shown in the fig. [6j the 
anisotropy of the Universe defined by (H — 7i)/(H + %) decreases like t . 



2 Semiclassical evolution in LQC 

2.1 Holonomy correction to Ashtekar variables 

As shown in the previous section, the massless vector field domination gives some predictions similar 
to the LQC. Nevertheless, it is worth to investigate the impact of vector fields on the LQC diagonal 
Bianchi I model - we shall compare it with results from GR and investigate more realistic models 
inspired by cosmic inflation. In this section we shall analyse the influence of the LQC holonomy 
correction on the classical evolution of scale factors, scalar field and a vector field. We shall focus 
on the effective, semi-classical equations of motion. Thus, we treat LQC as a modified, but still 
classical theory of gravity, alike f(R) theory. This means, that some effects, which come from a 
quantum nature of LQC, may be missing. We disregard the inverse triad correction, so only c\ and 
C2 are modified 

d ->■ —s'm(piCi) , (2.1) 
Pi 

where pi = Pi(pi,P2) will be specified later on. The effective gravitational Hamiltonian is now given 
by 

l 



T Ply/ Pi 



(2.2) 



2 i- \ Pi . o • /- \ ■ (- \ PlPZ 

(Pici) ^ + 2sm( / uici)sm {P2C2) 

p( P1P2 
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One shall remember, that the loop correction changes the physical interpretation of Cj. From now on 
ci ^ ja and C2 7^ 76. The low energy limit of the quantum theory shall give the General Relativity. 
Indeed, for the faci <C 1 one obtains 



1 

Pi 



sin(/ijCj) -> Ci 



C 



eff 



C 



grav • 



(2.3) 



The C ma t remains unchanged under the loop correction, since it is independent of Cj. The equation 



mat 



now takes the form C 



'eff 



+ C 



mat 



0. 



2.2 The old \i scheme 

In the isotropic LQC model one defines p, as y/A/p. To analyse the Bianchi I model one needs to 
introduce a generalisation of Ji. There are two models of LQC Bianchi I discussed in literature. The 
first one, described by e.g. [H] defines JJ14 as /2j = y/A/pi, where A is the eigenvalue of the area 
operator |24J. From the scalar constraint one obtains 



A7 2 



sin 



y/Acx \ p\ 



+ 2 sin 



Aci 



Pi / Pi 



sin 




(2.4) 



Equations ( 1.4|l.5 ) are unchanged under the loop correction and since p\ = ab(H+Ti) andp2 = 2a 2 H 
they give the Hubble parameters 



H 



A7 



sin 



H 



A7 



sin 



Aci \ 



Pi J 



cos 



Aci 



cos 




COS 




To obtain the Big Bounce we need H = % = at the Planck scale. There are two limiting cases, 
which are of interest: the low energy limit which we have already discussed, and the Planck energy 
limit with significant corrections from LQC. In the latter case it is easy to check, that 



sin 



Aci 



sm 



Ac 2 



1 



H = H = . 



(2.7) 



The critical energy density which corresponds to the Big Bounce is 



Pc 



1 



A7 2 I L\a cr 



+ 2 



Lob, 



20 CT 



(2i 



where a cr and b cr are the values of scale factors during the Big Bounce. One obtains the minimal 



value of the critical energy density for L\a cr = L2b cr , which is p c 



3 
A7 2 



The choice of the physical 



6 There are two values of p cr , which are often considered in literature: p cr ~ 0.82G 2 and p cr — 0.41G 2 . To find 
more about this issue see [TH]. For all numerical calculations we have chosen the p cr ~ 0.82G 2 , which means, that 
A ~ 3/( 7 2 0.82G 2 ) 
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length scales related to values of pi(t cr ) and P2(t C r) is arbitrary, for simplicity we set them to be 
equal. Unfortunately other choices give different evolution of pr£) and P2(t), which is a problem of 
anisotropic models of LQC in the old Jli scheme. The dependence of Bianchi I evolution on the 
choice of the fiducial cell and its rescaling is discussed in e.g. [TT| ITS]. 



Let us study effective equations of motion for c\ and 02- From eq. (1.5) one obtains 



1 



ci 



'Ac, . 
cos sin 



2A 7 



sin 



VP2 J Pi 



+ j(PA-p<p)VPl > ( 2 -9) 



(>2 



A7 




1.2 v^ci \ p\ . ( VAci \ [ v 7 AC2 
- sin — — r- -—77 + sin — ■ cos 



7(3PA+^) -(2.10) 



Those equations together with eq. ( 2.5|2.6 ) and equations of motion for matter fields give us the 
complete dynamics of the Universe. 



2.3 The new /j, scheme 

The problem of arbitrariness of the scaling in LQC of Bianchi I model may be solved by the redefi- 
nition of the /2j. Let us consider Ji\ = y 7 A/p2, p,2 = V^P2/pi [I3|T8]. Then one obtains p cr = -^1, 
so the critical energy density does not depend on the size of the fiducial cell. Form the Hamilton 
equations one finds 




(2.11) 
(2.12) 

(2.13) 

(2.14) 



One can see, that this choice of fii functions changes not only the dependence on the size of the 
fiducial cell, but also dynamics of Ashtekar variables. Nevertheless, as we will show, evolution 
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of matter field and scale factors are similar in both schemes. From the scalar constraint and eq. 



(2.11 2.12) one obtains 



n 



-H + 



3p 2 + 4H 2 p c 



3p 2 + 



2H 2 + p + p x 



3{H 2 -p) 



-11 

Per 



+ 4ff 2 + p 




3(H 2 -p) 

Per 
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3{H 2 -p) 



1 



Pc 



(2.15) 



which is the analytical relation between all measurable quantities of the model. It is worth to note, 
that for p ~ p cr one obtains T~L — —H, which is confirmed by numerical results. 

Let us discuss several scenarios for the vector field contribution to the evolution of the Universe. 
Usually one considers the isotropic Universe with subdominant vector field during inflation or ra- 
diation domination, or with the vector field domination (or its significant contribution) during the 
vector field oscillation phase. In following sections we shall analyse extensions of those model to the 
anisotropic Universe with the vector field domination and LQC corrections for both pi schemes. 



2.4 Massless vector field domination 

Referring to the LQC with a mass less s calar field domination, let us consider a massless vector field 
in the LQC frame. Since the eq. (1.15) is still valid one obtains p = p oc p^ 2 - Numerical solutions 
for the scale factors are shown in the fig. [TJ Ashtekar parameters in fig. [2j Hubble parameters 
during the Big Bounce in fig. [3] and the low energy limit for both Hubble parameters [3j The main 
difference with the classical scenario is the behaviour of b(t) near the Planck scale. The contraction 
phase starts from t = t cr , thus one obtains Kasner-like solution for t 3> t c 



2.5 Massive vector field domination 

The massless vector field domination leads to huge anisotropies and to contraction along the z 
axis. To consider a more realistic scenario let us assume the massive vector field domination. We 
have assumed kinetic initial conditions for the vector field evolution, so for p ~ p cr one obtains 
solutions similar to the massless scenario. On the other hand the growing value of A(t) stops the 
Kasner-like evolution, so b cr is only the local maximum of b(t). Numerical calculations show, that 
for m = W~ 5 m p i one obtains % > for t > 1000 or t > 1500 (old and new pn scheme respectively). 
The situation described in the subsection |1.3| is a good low energy limit for the LQC evolution, so 
the decrease of anisotropy shown in the fig. [6] remains valid in both schemes. 



3 LQC implications for realistic cosmological models 

So far we have studied the vector field domination in the early Universe. However, such scenarios 
can not produce inflation and lead to huge background anisotropies. Hence, let us turn to scenarios 
motivated by realistic inflationary models. 
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3.1 Massive vector with a massive scalar 



Let us consider the Universe dominated by a vector field Aa(t) with mass m and a scalar field 
(j)(t) with mass M. First of all we assume, that they are not coupled and that mass terms 
are time independent. Then one can introduce an a parameter defined by pAt = t cr ) = ap cr , 
PA(t = t cr ) = (1 — a)p cr . For numerical calculations we have assumed a = 0.1 (the vector field 
domination during the Big bounce) and a = 0.5 (vectors and scalars with the same contribution to 
Per)- I n both cases we assume m = M = 10 _5 m p ;. The evolution of the scale factor b(t) is shown 
in the fig. [4j The ratio pa/ Pcj> and evolution of % close to the Big Bounce are shown in fig. [5] 
The evolution of b(t) for both values of a shows, that the presence of the scalar field reduces the 
period of Kasner-like evolution. On the other hand the vector field dominates for the initial period 
t ~ O(10 3 t p j), which leads to significant background anisotropies during that period. 



3.2 Vector mass from the Higgs mechanism 

Let us consider the vector field with the time dependent mass term produced by the Higgs mechanism 
together with a constant mass term. This strong coupling is motivated by the desired evolution of 
the vector field in the inflationary regime (see eq. ( |1.11[ ) and discussion around it). The action of 
matter fields looks as follows: 

S = J d A x^~g {^R - + \d^4> + l 2 M) + m 2 )A^ - V((f>) \ . (3.1) 

It is easy to show, that the scalar field domination and g(4>) = — |V(^>) give the slow-roll solution 
for U(t) during inflation and the flat power spectrum of the initial inhomogeneities [I]. This makes 
it a useful choice, since often we would like the vector field to survive inflation and dominate the 
Universe after reheating. From the action (3.1) one obtains the following equations of motiorj^] 



4>+(2H + n)^> + V =0 , A + {2H -H)A + (m 2 -^V)A = . (3.2) 

At the level of Einstein equations the strong coupling term contributes to pa, so now = ^4> 2 — V 
and pa = — (m 2 — 2V/3)A 2 ). Let us focuse on the LQC corrections to the evolution of this 

model. Since both, vector and scalar fields, are effectively massless close to the Big Bounce, then 
the strong coupling shall not have a significant influence on their evolution in that energy range. 
The evolution of scale factors is shown in the fig. [6] The main deference appears for later times, 
when the extra mass term saves A(t) from domination. Coupling to the scalar field has been 



introduced to cancel time dependent part of the effective mass for the vector in ( 1.11 ) and the result 
is that the amplitude of the vector field does not decreases faster than that of the scalar field. This 
can lead to the strong backreaction of the vector field on the evolution of space-time, which would 
give large background anisotropies. To avoid that one needs U* < 10~ 3 , where U* is the value of 
the physical vector field at the moment of horizon crossing of its perturbations. Larger U+ would 
produce measurable vector perturbations in the CMB. 



7 Equations of motion can be also obtained from the scalar constraint ( 1.3 \ with m 2 — > m 2 + g((f>) 
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The other way to produce m? ~ —2H 2 is to consider a vector field with a non-minimal coupling 
to gravity. Unfortunately this model faces many problems in the LQC frame and so far we do not 
know how to to obtain effective equations of motion in this case. 



4 Conclusion 

In the initial section of this paper the evolution of an anisotropic Universe dominated by vector field 
in GR frame has been studied. We have obtained the analytical solution for the massless vector 
field domination. In particular, it has been found, that there appears the quasi-Big Bounce solution, 
which is characterized by a finite value of the maximal energy density (p(ti) = pi < oo), lack of 
initial curvature sing ularity (R = and W v R^,W uaP R^p < oo) and a non-zero initial value 
of the scale factor a (a(if) = ai ^ 0) or H{ti) = 0. On the other hand, along the direction set 
by the vector field one obtains characteristic features of the Big Bang for t — >• ij. The maximal 
energy density pj is a constant of integration. The b(t) grows for t £ (tj,tj + 8/(3^fpj)), and after 
that period one obtains Kasner-like solution. For the massive vector field domination one finds the 
isotropic Universe at late times. 



In the section [2] the LQC holonomy corrections to the simplified Bianchi I with a combination 
of vector and scalar fields have been considered. The effective equations of motion for Ashtekar 
variables, see sec. 2.1 2.3 have been found for both, old and new pi schemes. We have performed 
numerical simulations of the evolution of matter fields and scale factors for following scenarios: (i) 
massless vector field domination |2.4| we have obtained the Big Bounce at t = t cr and Kasner-like 
solution with decreasing value of b(t). Thus b(t cr ) is the global maximum of b(t). (ii) Massive vector 
field domination, see 2.5 one obtains Kasner-like solution for the first t ~ O(10 3 t p ;), later on the 
Universe expands in all directions. 



In the final part of the paper we have discussed scenarios inspired by realistic inflationary models. 



Domination of a combination of a vector and a scalar field, see 3.1 produces the transition between 
contraction and expansion along the z axis. Regardless of the initial contribution to the energy 
density, the vector field dominates the Universe, but one obtains the scalar field domination long 
before the vector oscillation phase. The relatively strong coupling between the massive vector and 
the massive scalar gives b > for t ~ O(0. lt p [). It is interesting to note, that non-trivial features 
on the evolution of the Universe appear in this case for energies much lower than the Planck scale. 
Also, it is important to note, that even during the b < phase our model is consistent with the 
LQC, since both Ashtekar parameters reach their global minima at t cr . One shall also note, that 
both forms of pi functions gives similar results even at high energies, but this similarity can be 
at most qualitative, because of the scaling problem in the old p scheme. Future analysis shall be 
focused on the birth of quantum fluctuations in the anisotropic Universe with the LQC holonomy 
correction and on the non-minimal coupling of a vector field to gravity. Further research is needed 
to find out, how generic the results are. 
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Figure 1: In the panel (a) we present the evolution of a(w)/ai for massless vector field domination 
in GR. Panel (6) shows the evolution of b(w)/{aiE). The red line corresponds to b(w), and the blue 
dashed line corresponds to b(w) oc u) -1 / 3 , which is the low energy limit for b(w) evolution. Panels 
(c) and (d) show a{t)/a cr and b(t)/b cr as a function of time in Planck units for the LQC with the 
massless vector field domination. Blue and red dashed lines corresponds to the old and the new jli 
scheme respectively. One can see, that b cr is the global maximum of b{t). 
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t 



Figure 2: Left, middle and right panels show pi, p2 and V respectively for the Universe dominated 
by the massless vector field. Dashed red and blue lines denotes new and old fii scheme respectively. 
It is worth to note, that even if 6(t cr ) — b max Ashtekar variables p\ and P2 have their global minima 

at t = try = 0. 




Figure 3: All panels show H and % for the LQC with a massless vector field domination for the old 
(blue line) and the new (red dashed line) Jli scheme. For t > one founds % < which means, that 
b cr is the global maximum of b(t). Panels (c) and (d) show the LQC solution (blue and red lines) in 
a wider time range together with the GR limit (dashed yellow lines). In the panels (c) and (d) the 
moment of the Big Bounce is not shown for the LQC curves due to assumed scale of the plot. 
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Figure 4: Left and right panels show b(t) for LQC with the massive scalar and vector fields for 
a = 0.5 and a = O.lrespectively. Blue and dashed red lines denotes old and new fa scheme. Initially 
b(t) decrees, but after a while the scalar field starts to expand Universe in that direction. The a(t) 
grows continuously from Qj — O'er 

like in the vector field domination scenario. One obtains very 
similar results for the vector mass coming from the Higgs mechanism. 




Figure 5: Left panel shows pa/ P4> f° r LQC with the massive scalar and vector fields for a = 0.1 (blue 
and red lines) and a = 0.5 (yellow and green lines). The old and the new fa scheme are denoted 
as continuous and dashed lines respectively. The vector field initially dominates the Universe, but 
after 10 3 Planck times for a = 0.5 and after 4 x 10 3 Planck times for a = 0.1 one obtains scalar 
field domination. The right panel shows the evolution of %. For a = 0.1 one obtains % > when 
t > 0.5t p i 
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Figure 6: In the right panel we present the evolution of H (blue and red lines) and % (yellow 
and green lines) for domination of a scalar and a vector field with a strong coupling in the LQC 
frame. The old and the new /2j scheme are denoted as dashed and continuous lines respectively. 
For the numerical calculations we have assumed purely kinetic initial conditions. In the left panel 
we present the evolution of the Universe anisotropy during the period of the domination of the 
oscillating vector field. Dashed line denotes the analytical solution for (H — T~L)/{H + T~L) for the 
pressureless anisotropic Universe. This solution corresponds to p ^ p cr , where both schemes of fli 
gives the same results 
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